Manipulating Majorana fermions in one-dimensional spin-orbit coupled atomic Fermi 

gases 
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Majorana fermions are promising candidates for storing and processing information in topological 
quantum computation. The ability to control such individual information carriers in trapped ultra- 
cold atomic Fermi gases is a novel theme in quantum information science. However, fermionic atoms 
are neutral and thus are difficult to manipulate. Here, we theoretically investigate the control of 
emergent Majorana fermions in one-dimensional spin-orbit coupled atomic Fermi gases. We discuss 
(i) how to move Majorana fermions by increasing or decreasing an effective Zeeman field, which acts 
like a solid state control voltage gate; and (ii) how to create a pair of Majorana fermions by adding 
a magnetic impurity potential. We discuss the experimental realization of our control scheme in an 
ultracold Fermi gas of 40 K atoms. 

PACS numbers: 03.75.Ss, 71.10.Pm, 03.65.Vf, 03.67.Lx 



Majorana fermions - particles that are their own an- 
tiparticles - were proposed by Ettore Majorana in 1937 
[1(. They are finding ever wider relevance in modern 
physics [2|. A well-known proposal is the use of Majo- 
rana fermions as information carriers in topological quan- 
tum information processing [3, 4]. Due to the robustness 
of topological quantum computation, research into the 
creation and manipulation of Majorana fermions has be- 
come an active research topic in a variety of fields of 
physics, ranging from condensed matter physics to ultra- 
cold atomic physics. Most theoretical schemes of process- 
ing topological quantum information based on Majorana 
fermions have relied on using solid-state systems [B|, @ . 
Evidence of the existence of Majorana fermions in hy- 
brid superconductor-semiconductor nanowire devices has 
been reported very recently Q- 

The realization of Majorana fermions in ultracold 
atomic systems has also been discussed extensively 0- 
13| . Compared with solid state systems, a unique ad- 
vantage of ultracold atoms is their unprecedented con- 
trollability and purity, which may be expected to lead 
to reduced decoherence. A vast range of interactions, 
geometries and dimensions are possible: using the tool 
of Feshbach resonances, one can control the interatomic 
interactions very accurately. By using the technique of 
optical lattices, one can create artificial one-dimensional 
(ID) or two-dimensional (2D) environments to explore 
how physics changes with dimensionality. In this con- 
text, topological quantum computation based on Majo- 
rana zero-energy modes in the vortex core of a 2D atomic 
Fermi gas has been discussed @. However, it is very diffi- 
cult to experimentally generate and manipulate a vortex 
lattice in a fermionic superfluid. 

In this Brief Report, we theoretically investigate the 
manipulation of Majorana fermions in one-dimensional 
spin-orbit coupled atomic Fermi gases. The existence of 



Majorana fermions in such systems was recently checked 
by us, by performing fully microscopic calculations within 
the mean-field Bogoliubov-de Gennes (BdG) theory |{|. 
Here, we study how to control Majorana fermions using 
a background magnetic field combined with a magnetic 
impurity potential. Our investigation is motived by the 
recent realization of a three-dimensional spin-orbit cou- 
pled Fermi gas of 40 K and 6 Li atoms, at ShanXi Univer- 
sity 14 1 and at MIT , respectively. We anticipate that 
a ID spin-orbit coupled Fermi gas will be created very 
soon, by using 2D optical lattices |16j. At the end of this 



paper, we discuss the possible experimental realization of 
our control scheme for 40 K atoms. 

Model Hamiltonian. - Our theoretical approach is 
outlined in detail in the previous work Here we ex- 
tend this approach to include a classical magnetic impu- 
rity. In the following, we introduce briefly the essential 
ingredients of the theory. The model Hamiltonian of a 
ID atomic Fermi gas with spin-orbit coupling, as realized 
at ShanXi Q and at MIT [lj| , can be written as 

U = j dx^ (x) [U s a (x) - ha z + Xka y ] tp (x) 

+»*/**?(.)*!(.)*(.)*(.). (i) 

where ip^ (x) = [ipj. (x) , -0j {x)] denote collectively the 
creation field operators for spin-up and spin-down atoms, 
h denotes the strength of an effective magnetic field caus- 
ing a Zeeman splitting, while Xka y = —i\(d/dx)a y is the 
spin-orbit coupling term with coupling strength A. In ad- 
dition, g v and a z are the usual 2x2 Pauli matrices while 
gio = —2h 2 /{maio) is the effective Hamiltonian cou- 
pling given an s-wave scattering length in one dimension 
of am- We note that the existence of spin-orbit coupling, 
magnetic fields and quantum transport to the cloud edge 
remove all of the usual symmetry properties of Fermi 
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gases, and these systems therefore provide an opportu- 
nity to investigate the least symmetric of the nonstan- 
dard symmetry classes in normal-superconducting hybrid 
structures (l7l|. 

The single-particle Hamiltonian H^(x) describes 
the atomic motion in a harmonic trapping potential 
muj 2 x 2 /2, together with a scattering potential of a clas- 
sical magnetic impurity, modeled as — Vi mp (x) a z . Com- 
bining these terms, %^{x) takes the form, 

H%(x) = -(h 2 /2m)(d 2 /dx 2 )+(m/2)cj 2 x 2 -V imp (x)a z -», 

(2) 

where fi is the chemical potential. The impurity scatter- 
ing potential with strength Vq is taken to have a gaussian 



shape, so that 

V imp (x) = [Vo/(V2^d)} exp[-(x - x ) 2 /(2d% (3) 

where xq and d are the position and width of the impurity 
potential, respectively. Thus, in the limit of infinitely 
narrow width where d — > 0, we obtain a Dirac delta- 
function like potential, Vi mp (x) = Vq5{x — x$). 

We solve the model Hamiltonian ([1]) using a 
mean-field Bogoliubov-de Gennes (BdG) approach. 
The wave-function of low-energy quasiparticles 

[u tv (x) , u iv (x) , v tv (x) , v iv (x)] T with energy E v 
satisfies the BdG equation, 
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where the pairing field is defined as: 

A(x)=-( 9lc /2)^K(x)^(x)/(S,)+ 
v 

+ u u {x)v^(x)f{-E v )} (5) 

and the density of spin-a atoms is given by n a (x) = 
(1/2)E,[I^(.t)| 2 f(E n ) + \v ari (x)\ 2 f(-E n )}. Here, 
f(E) = l/[e B /( feBT ) + 1] is the Fermi distribution func- 
tion at temperature T. The BdG equation Q can be 
solved self-consistently using a hybrid method detailed 
in Ref. [3. 

In the following calculations we consider a typ- 
ical dimensionless interaction parameter of 7 e 
(aho/aw) I (7T./V 1 / 2 ) = f .6, where a^o = \fh/(mw) is the 
characteristic oscillator length of the trap and N = 100 is 
the total number of fermionic atoms. For spin-orbit cou- 
pling, we take Xkp/Ep = 1, where kp and Ep are respec- 
tively the Fermi wave-vector and Fermi energy. The unit 
of length is the Thomas-Fermi radius, x F = N^ 2 a ho . All 
calculations will be carried out at zero temperature, as 
the inclusion of a finite but small temperature will only 
marginally affect our results. 

A brief summary of previous results. - In the 
absence of magnetic impurity, the spin-orbit coupled ID 
Fermi gas in a harmonic trap will locally enter an in- 
teresting topological phase when the Zeeman field satis- 
fies h > h c (x), as we reported earlier in Ref. @. Here 
h c (x) = \J \l 2 (x) + A 2 (x) with \x (x) = fj, — muj 2 x 2 /2 
is the local critical field at position x. For interaction 
strength 7 = 1.6 and spin-orbit coupling Xkp/Ep = 1 
mentioned above, this occurs at h > h c ~ 0.35Ep See, 



for example, the phase diagram in Fig. 1. At higher 
magnetic fields of h > h ci a mixed phase emerges, con- 
sisting of a standard BCS superfluid at the trap center 
and a topological superfluid at the two edges of the trap 
(see Fig. 2). 

As a salient feature of topological superfluids, Majo- 
rana fermions appear at the edges. These Majorana 
fermions are nothing but the zero-energy solutions of the 
BdG equation. Because of the particle-hole symmetry 
of the BdG equation, the field operator of Bogoliubov 
quasiparticles satisfies = T 1 _ E . If the energy E = 0, 
we would then have Tq = Tj, which is the defining fea- 
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Figure 1: (color online) Phase diagram determined from the 
lowest energy of Bogoliubov quasiparticle spectrum. The 
inset shows the contour plot of local density of state at 
h/Ep — 0.5. The contributions from Majorana fermions are 
highlighted by circles. 
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ture of Majorana fermions. Related to this zero-energy 
feature, it is easy to see that spatially resolved radio- 
frequency (r.f.) spectroscopy, which measures the lo- 
cal density of state of the atoms, defined as p(x, E) = 
(1/2) Y Jav [\u.r,{x)\ l 5{E-E ri ) + \v (7ri {xf S(E+E ? )], pro- 
vides a useful means to identify Majorana fermions. In 
the inset of Fig. 1, we can identify from r.f. spectroscopy 
four zero-energy Majorana fermions at the Zeeman field 
h = 0.5Ep, two for each topological superfluid. 

Manipulating Majorana fermions. - Let us con- 
sider how to manipulate Majorana fermions in ID atomic 
Fermi gas. We focus on two operations: to move and to 
create Majorana fermions. 
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Figure 2: (color online) Positions of Majorana fermions as a 
function of Zeeman field. 

We propose to use the Zeeman field, as an analog to a 
voltage gate in solid-state systems 0|, to transport Ma- 
jorana fermions through the whole Fermi cloud. Indeed, 
the position of Majorana fermions depends critically on 
the Zeeman field, as shown in Fig. 2, where we plot the 
boundary of each superfluid phase. The four Majorana 
fermions sit exactly on these boundaries. Therefore, by 
adiabatically increasing or decreasing the Zeeman field, 
the inner two Majorana fermions will move towards or 
away from the trap center. The position of the outer 
two Majorana fermions, however, is less affected by the 
Zeeman field. 

To create new Majorana fermions in the system, we 
may consider either (i) creating a new topological super- 
fluid in some areas of the Fermi cloud or (ii) introducing 
a strong defect in the existing topological superfluids. In 
the former case, a new pair of Majorana fermions will 
appear at the edges of the induced new topological su- 
perfluid. In the latter case, the strong defect will create 
a low-density hole, which is basically vacuum and thus 
topologically trivial. New Majorana fermions will then 
be trapped at the hole center. A well-known example of 
the latter case is a vortex in 2D topological superfluid, 
which hosts a zero-energy Majorana fermion mode inside 
the vortex core. 

The main result of this work is that these two ideas (i) 



and (ii) of creating new Majorana fermions can be real- 
ized by introducing a classical magnetic impurity. The 
magnetic impurity serves as an attractive or repulsive 
scattering potential for spin-up or spin-down atoms, re- 
spectively, as we mentioned earlier in the model Hamil- 
tonian. Due to the finite width of the impurity scatter- 
ing potential, the first idea (i) seems more practical. In 
the following, we will focus on this idea and examine it 
by using self-consistent BdG calculations. To test this 
concept, we will put an impurity at the center of the 
Fermi cloud (i.e., xq = 0) and take a potential width 
d = Q.OTxf- The scattering strength of the impurity is 
taken as Vq = O.ISxfEf . 

The existence of impurity induced Majorana fermions 
is clearly evident from the local density of state, which 
is measurable using spatially resolved r.f. spectroscopy, 
as reported in Fig. 3a using a linear contour plot. For 
comparison, we show also in Fig. 3b the expected spec- 
troscopy scan without an impurity. It is clear that 
a new pair of Majorana fermions is created at about 
x = ±0.lxp, while other four Majorana fermions almost 
keep the same positions as those without the impurity po- 
tential. The creation of new Majorana fermions can be 
easily understood. The magnetic impurity potential with 
a relatively broader potential width may be regarded as 
a local Zeeman field. By tuning the strength of the mag- 
netic impurity, we can greatly enhance the local Zeeman 
field at a particular area. This creates a local topolog- 
ical superfluid that hosts new Majorana fermions at its 
edges. To check this picture, we have calculated the cri- 
terion for a local topological superfluid, h — h c (x) > 0. 
As highlighted in Fig. 4a by a pink cross-pattern, we find 
that this criterion is indeed satisfied at about the origin, 
where the magnetic impurity is located. The superfluid 
gap is significantly reduced at the magnetic impurity. 

We finally check how the density distribution is af- 
fected by magnetic impurity. In Fig. 5a, we present 
the spin-up and spin-down density distributions and their 
difference An(x) = n^(x) — ni(x). Compared with the 
result without the impurity (Fig. 5b), we find that the 




Figure 3: (color online) Linear contour plot of the local den- 
sity of states with (a) or without (b) a magnetic impurity at 
h/E F = 0.5. 
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Figure 4: (color online) Criterion for topological superfluid 
h — h c (x) (solid lines) and the superfluid order parameter 
A(i) (dot-dashed lines) at h/Ep — 0.5. The upper and lower 
panels show the results with and without a magnetic impu- 
rity, respectively. The cross-patterns highlight the areas of 
topological superfluids. 
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Or, which flip atoms from \l) to |t) spin states and si- 
multaneously impart momentum 2hkn (or vice versa), 
via the two-photon Raman process [l4j • This is the same 
technique used in the National Institute of Standards and 
Technology (NIST) to generate synthetic spin-orbit cou- 
pling in a Bose- Einstein condensate of 87 Rb atoms (l9| . 

In our model Hamiltonian, the spin-orbit coupling 
strength is related to the wave-vector k^ by A = 
h 2 kn/ (2m), and the effective Zeeman field can be tuned 
by the strength of Raman laser beam, i.e., h = Qr/2. 
To make the system one-dimensional, one can impose 2D 
optical lattices that restrict atoms moving along a single 
tube. This was already demonstrated at Rice University 
to obtain an unbalanced ID Fermi gas [lij]. It is possi- 
ble to create an impurity potential by using a laser beam 
that has a sufficient narrow beam width, as schemati- 
cally shown in Fig. 6. By suitably tuning its frequency, 
the scattering potential caused by the laser beam can be 
attractive for the spin-up state and repulsive for the spin- 
down state. This localized probe therefore behaves like a 
classical magnetic impurity. All the techniques required 
to simulate the model Hamiltonian Eq. ([1} are therefore 
within the reach of current experiments. 

We note that our scheme does not yet lead to a fully 
developed quantum gate, which in any event requires a 
careful analysis of quantum dynamics during the switch- 
ing process. Instead, it provides a simple demonstration 
system whereby the concepts of creating and moving Ma- 
jorana fermions in a quasi one-dimensional environment 
- like 'birds on a wire' - can be usefully investigated. 

Conclusions. - In summary, we have theoretically 
proposed how to manipulate Majorana fermions in a one- 
dimensional spin-orbit coupled atomic Fermi gas in a har- 
monic trap. The transportation of Majorana fermions 
through the Fermi cloud may be controlled by an effec- 
tive Zeeman field, while the creation of pairs of Majorana 
fermions can be achieved through an artificial magnetic 
impurity created using a dimple laser beam. Our pro- 



Figure 5: (color online) The spin-up and spin-down density 
distribution, (dashed lines) and (solid lines), and their 
difference An = rif — nj. (dot-dashed lines) at Ji/Ef = 0.5. 



local density is enhanced and reduced for the spin-up and 
spin-down atoms, respectively, due to the magnetic im- 
purity. This is consistent with the picture of inducing a 
local Zeeman field, as mentioned in the above. 

Experimental proposal. - The model Hamiltonian 
Eq. ((T|) can be realized using ultracold 40 K atoms 
The two spin- 1/2 states can be chosen as the two mag- 
netic sublevels |f) = \F = 9/2, m F = -9/2) and \i) = 
\F = 9/2, m F = -7/2) of the F = 9/2 hyperfine level. 
To create the spin-orbit coupling, one can use a pair of 
counter-propagating Raman laser beams with strength 
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Figure 6: (color online) Schematic plot of the experimental 
setup. The spin-orbit coupling is induced by a pair of Raman 
laser beams, while the magnetic impurity with finite width is 
created by a dimple laser beam. 
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posal may be viewed as the cold-atom analog of the con- 
trol gates proposed in the solid-state systems Q. We 
anticipate that these control operations can be realized 
using ultracold 40 K atoms Q. 

We thank Hui Hu and Han Pu for useful discussions. 
This work is supported by the ARC Discovery Project 
(Grant No. DP0984637) and the NFRP-China (Grant 
No. 2011CB921502). 
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